Selection effects on GRB spectral-energy correlations 
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Abstract. Instrumental selection effects can act upon the estimates of the peak energy E°^ k , the fluence F and the peak flux P 
of GRBs. If this were the case, then the correlations involving the corresponding rest frame quantities (i.e. £p ea k, £; so and the 
peak luminosity L; so ) would be questioned. We estimated, as a function of Ep^, the minimum peak flux necessary to trigger a 
GRB and the minimum fluence a burst must have to determine the value of E°^ k by considering different instruments (BATSE, 
Swift, BeppoSAX). We find that the latter dominates over the former. We then study the £°k| k -fluence (and flux) correlation in 
the observer plane. GRBs with redshift show well defined E^^-F and Ep^-P correlations: in this planes the selection effects 
are present, but do not determine the found correlations. This is not true for Swift GRBs with redshift, for which the spectral 
analysis threshold does affect their distribution in the observer planes. Extending the sample to GRBs without z, we still find 
a significant E^^,-F correlation, although with a larger scatter than that defined by GRBs with redshift. We find that 6% are 

outliers of the Amati correlation. The E^^-P correlation of GRBs with or without redshift is the same and no outlier is found 
among bursts without redshift. 
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INTRODUCTION 

Thanks to the discovery of the afterglow emission in long Gamma-Ray Bursts (GRBs), redshift measurements of these 
sources became available. The knowledge of GRB distances made possible to estimate their intrinsic properties and 
several correlations between prompt and afterglow properties were found. One of the most intriguing is the so called 
Ghirlanda correlation (111]), linking the rest frame peak energy £p ea k to the collimated energy Ey emitted during the 
prompt. The possibility to use this correlation for cosmological purposes ([2]) makes it very appealing. However, the 
possible presence of selection effects, the lack of a physical interpretation, the difficulty in calibrating the correlation 
and other problems mainly related to the jet break time identification, raised objections about the validity of this 
correlation and/or its application in cosmology. The lack of outliers (except for GRB980425 and GRB031203, but see 
JH) and the small dispersion of points around the best fit line are encouraging, but at this stage the paucity of GRBs 
defining this correlation (27 GRBs) prevents us from drawing any firm conclusion. 

The paucity of GRBs with known £p ea k and Ey is partially due to the difficulty to estimate the jet break times. 
They are expected in the afterglow light curves at late times (about 1 day), when typically the light curve is not well 
sampled. At late times, the observations are lacking and difficult (e.g. for the presence of the host galaxy and/or a 
possible contamination from the underlying supernova emission) and they are not sufficient to reveal the presence of a 
break in the powerlaw decay behavior. The existence of other breaks and, more generally, the very complex temporal 
behavior discovered both in the X-ray and optical light curves, makes the situation even more confused. The different 
temporal decay of the X-ray and optical light curves raised the question of a different origin for these two emissions. 
101 suggested that often the standard external shock model only accounts for the optical emission, whereas the X-ray 
has another (internal) origin. According to this interpretation, Jj] caution about the identification of jet break times, 
pointing out that they should be identified in the optical light curves, at late times. According to the modeling of @], 
these breaks could also be chromatic. 

When the jet break time is not available, we can only estimate the isotropic equivalent energy £i S0 , i.e. the energy 
emitted during the prompt by assuming an isotropic emission geometry. The quantity £i so still shows a correlation with 
£peak (the so called Amati correlation, J7[]), in the form fipeak^^iso ' 5 ■ This correlation, found before the Ghirlanda one, 
displays a larger scatter and a different slope. Assuming the validity of the E pea k-Ey correlation, the presence of this 
large scatter in the Amati correlation can be easily explained, assuming that a wide range of jet opening angles 9j can 
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FIGURE 1. £peak"£iso an d ^peak _ Aso correlations for 83 GRBs with measured redshift and spectral parameters. The slopes 
obtained from the fit with a power law model are respectively 0.48 ±0.03 and 0.40 ±0.03. Modeling the scatter distributions with 
a gaussian function we found respectively a = 0.23 and a = 0.28. 



correspond to bursts with same Ey (and therefore, for the E pel ±-Ey correlation, same £p ea k)- This range of 9j values 
reflects into a range of Ei so values, for the same given £p ea k- 

Independently from its origin, the large dispersion of the Amati correlation makes it unsuitable for the cosmological 
applications, unless systematics, unknown, extra scatter terms are introduced. 

Despite these problems, the Amati correlation appears robust. Although the number of GRBs with known £p e ak and 
Ei so is reasonably large (83 objects inja]), only two of them do not obey the Amati correlation. Note that these two 
outliers are very peculiar GRBs (see yfl). Moreover, the validity of the correlation has been extended to XRFs and 
now spans five orders of magnitudes in Ei so , from XRFs to very bright GRBs. This behavior calls for a theoretical 
explanation. Some works (0], lloll . iflllo suggest that its existence is due to the presence of selection effects. To check 
this possibility it is necessary to perform a systematic and accurate study of all the possible selection effects, to better 
understand if and how they affect the correlation. 

Another correlation involving similar quantities is between £p ea k and L; so (Yonetoku correlation). fl2ll found this 
correlation defining Lj so as the luminosity emitted at the peak of the light curve. In its original form, this correlation 
appeared slightly tighter than the Amati one, but with a similar slope. We will also study selection effects on this 
correlation. Fig.[TJshows both the £p e ak-Eiso an d the E pea ^-Li so correlations updated to April 2008 ([8]). The analysis 
performed on this sample (83 GRBs) shows that the £p e ak-£iso correlation displays a larger scatter than the £p e ak-£iso 
correlation and therefore the latter seems to be more fundamental (see (8J for details). 
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FIGURE 2. Distributions of 83 GRBs with known redshifts in the £°k| k -Fluence and ^pe^-Peak flux planes. In both cases, they 
define a correlation in the observational planes. Shaded regions in the right panel show the ST (see text) for different instruments: 
for each instruments it is possible to perform a reliable spectral analysis only for bursts that lie on the right side of these regions. In 
both the panels, the curves show the TT (see text) for different instruments: only bursts at the right side can be detected. 



INSTRUMENTAL SELECTION EFFECTS 

We are interested in studying the rest frame correlations £peak~£iso ind £peak~^4so- However, instrumental selection 
effects act on the corresponding observational quantities: E°^ k , Fluence (F) and Peak-Flux (P). In Fig.[2]we report the 

sample of 83 GRBs with known redshift plotted in the observational planes E^^-F and E^^-P. Also in these planes, 
this sample of GRBs defines a correlation. 

Instrumental limits are obviously present. For example, very faint bursts can not be detected by a given instrument 
if they are below its trigger threshold. It is possible that the correlation is completely due to selection effects that make 
accessible only a small part (a stripe) of the whole observational plane. However, note that the region at the right side 
of the £°g a s k -F (^peajj-P) correlation can not be affected by instrumental selection effects, since bursts in this region 
would have an intermediate/low £p ea k (easily measurable by instruments) and a high fluence (peak flux). The lack 
of bursts in this region strongly suggests that bursts with these characteristics do not exist or they are very few. On 
the left side of the two correlations we can face with two different situations: the unavoidable selection effects acting 
in these regions can i) lie far from the position of observed GRBs (in this case the observed distribution of points is 
not determined by selection effects) or ii) lie very close to points, indicating that they determine the shape of GRBs 
distribution in this part of plane. To discriminate between these two cases, we need to quantify the selection effects. 

We identified the two most probable instrumental selection effects that can affect the sample: 

• the trigger threshold (TT): it is very complex to quantify what characteristic a burst must have to be detected by 
a given instrument. To a first approximation, we can affirm that it must have a minimum photon peak flux. This 
translate into a minimum energy flux, whose value depends on the spectral properties of the burst (especially 
Epeak)- F° r several instrument [ 13] derive the curves of minimum peak flux as a function of E°^ k . We plot these 
curves in the £°^ k -P plane, after accounting for the bolometric correction. Considering a typical conversion factor 
between fluence and peak flux of ~ 6 ( Il4ll ). we plot these curves also in the E^^-Fluence plane. 



• the spectral analysis threshold (ST): we need a minimum fluence to perform a reliable spectral analysis and 
determine E°^ k and the spectral shape. Hill consider this problem by accounting for the detector response and 
the typical background of each instrument. For each £°g a s k they determine this limiting fluence and obtain a curve 

in the E^^-F plane. Its position depends on the burst duration. In Fig.|2](left panel) we show the ST as a shaded 
region: the left (right) edge represents the curve obtained for bursts lasting 5 (20) seconds. 

Comparing the position of points in Fig. [2] detected by a given instrument with the corresponding threshold curves 
we can draw several conclusions. In the E°^l k -F plane the ST represents a more important selection effect compared to 
the TT. Swift bursts (filled squares) are strongly affected by the corresponding ST limiting curves. Note that for Swift 
bursts we mean bursts whose E°^ k has been determined from the modeling of the BAT data. However, this instrument 

can only measure E°^ k within the range 15-150 keV. Consequently, the Swift sample is distributed in a very narrow 

range in E°^ k , smaller than the scatter of the £°g a s k -F correlation as defined by all the pre-Swift bursts. For this reason, 
we caution about the use of Swift bursts to test the existence of this correlation. 

Our results also show that the pre-Swift GRB sample, containing a fraction of bursts detected by BeppoSAX and 
BATSE, is not affected by the corresponding limiting curves. This is particularly evident for high £°g a s k , where the 

limiting curves lie far from the points. A similar conclusion can be drawn for the ^p^-P correlation (right panel). 
Points are slightly affected by the corresponding TT in the region at low peak flux, but the behavior according to 
which higher peak fluxes correspond to higher E°^ k is not due to the threshold. 

In both planes however, the lack of bursts (with known redshift) with intermediate/large E°^ k and small fluence or 
peak flux (i.e. between the present sample of GRBs and the limiting curves) could be due to the additional presence of 
a still-not-understood selection effect, for example linked to the redshift determination. The first step to investigate this 
possibility is to search if there exist GRBs which populate the region on the left-hand side of the E°^-F and £pg ak -P 
correlations. This issue can be solved by including, in both plots, all the bursts with measured spectral properties, 
without the request of knowing their redshifts. 



GRBS WITHOUT REDSHIFT MEASUREMENT 

In this section we study the very same issue tackled in the previous one, but considering all bursts found in literature 
with measured £°g ak and spectral properties. For this reasons, we considered Hete-II bursts ( fl5ll ). Swift bursts (dgj]) 
and the sample of bright BATSE bursts (| 16]). Moreover, we collected the preliminary results arising from the analysis 
of the Konus-Wind spectra as reported in the GCN circulars (29 GRBs). To extend the sample of BATSE bursts to 
lower fluences we performed the spectral analysis of 100 BATSE bursts down to 10~ 6 erg/cm 2 (see [8] for details). 

Fig.[3](left panel) reports the distribution of all these bursts in the E°^ k -F plane. Different symbols refers to different 
instruments. We note that the entire sample of GRBs without redshift seems to occupy a larger region with respect to 
GRBs with redshift (filled squares) studied in the previous section. In particular, they enlarge the distribution toward 
the ST curves. However, moving towards lower fluence values their density decreases (see ^ for more details). In 
other words, points do not uniformly fill all the accessible plane, but are concentrated along a strip whose slope and 
position is not determined by the considered selection effects. Therefore, the existence of the Amati correlation must 
have a physical origin. However, the dispersion of the whole sample of GRBs is larger then that defined only by 
GRBs with redshift. This suggests that, if we were able to measure the distance of all GRBs in Fig. [3] we should 
found an Amati correlation with different properties, such as a larger scatter and, likely, a different slope. The grey 
region in the top left corner of the E^j^-F plane (see Fig. [3} shows the region where bursts are outliers of the Amati 
correlation (as defined by the sample of 83 GRBs with redshift). This means that GRBs that lie in this region can not 
be consistent with the correlation for any redshift. For 'consistent' we mean the burst must fall within the 3(7 scatter 
of the correlation. We found that the 6% of the sample is not consistent with the Amati correlation. 

The right panel of Fig. prefers to the ZSp^-Peak flux plane. Also in this case it is evident the presence of a 
correlation, i.e. of a distribution of points that preferentially lie along a strip. This behavior cannot be explained in terms 
of selection effects, since the shape and the position of the TT curves cannot account for the observed distribution of 
points. The difference with respect to the case of the E^^-F correlation is that, here, the distribution of points without 
redshift seems to better agree with that of GRBs with redshift. This suggests that in the future, when a larger sample 
of GRBs with redshift will be available, the properties of the Yonetoku correlation (slope and scatter) will resemble 
those of the present correlation. Also in this plane we have estimated the regions of outliers and we found only one 
burst which cannot be consistent with the correlation. 
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FIGURE 3. Left panel: bursts with (filled squares) and without (other symbols) known redshift in the Ep^-Fluence plane. 
Shaded curves represent the ST curves for BATSE and Swift/BAT instruments (see Fig. [2}. Shaded regions in the top left corner 
show which bursts can not be consistent with the 3o" scatter of the Amati relation defined in the rest frame plane by GRBs with 
redshift. Right panel: E°^-Peak flux plane. Curves show for several instrument the TT, i.e. the minimum peak flux required to 
detect the burst. In this case only one burst fall in the region of outliers. 



CONCLUSIONS 

To study the role of possible instrumental selection effects on the Amati and Yonetoku relation, we have considered the 
distribution of GRBs with and without known redshift in the observational £°^ k -Fluence and E°^ k -Peak flux planes. 
Following the analysis performed by 11411 . we refer to two different instrumental biases: the trigger threshold (TT, the 
minimum peak flux required to trigger a burst) and the spectral analysis threshold (ST, the minimum fluence needed 
to constrain the GRB spectral properties). These curves depend on E°^ k and define what part of the observational 
planes is accessible: only bursts at the right side of both the curves can be plotted in these planes, the other bursts 
have no sufficient flux to be detected or no sufficient fluence to recover £°g a s k and the spectral shape from the spectral 
analysis. We note that in both the observational planes GRBs are not spread in the region free from instrumental 
selection effects, but define a correlation. Note that the shape of this concentration of points is not determined by the 
considered thresholds. Their only effect is to cut the part of the correlations corresponding to the smallest E°^ k and 
Fuences/Peak fluxes. From the comparison between bursts with and without redshift we can conclude that there exists 
an £peak-£iso correlation not determined by selection effects, even if its slope and scatter may be different from what 
we know now. On the other hand, we suggest that, once a large number of bursts with redshift will be available, the 
Yonetoku correlation will preserve more or less the present scatter and slope. 

Another hint in favor of the relevance of the £p e ak-Aso correlation comes from short bursts. A detailed analysis of a 
large sample of short bursts performed by lfl7ll shows that they have a similar E°^ and Peak flux of long GRBs and, 

indeed, they populate the same region in the £ , °g a s k -Peak flux plane. This suggests they can be consistent with the same 
£peak-£iso correlation defined by long GRBs (this is also confirmed by the few short bursts for which we know the 
redshift). Short GRBs, instead, are inconsistent with the distribution of long GRBs in the E^^-Fluence plane, since, 

for the very same E°^ k they have a smaller fluence. This implies that the majority of short GRBs are outliers of the 
Zspeak-^iso correlation defined by long bursts. 
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